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Abstract 

We consider the 0{a'^) hadronic corrections to the energy spectrum of the decay 
electron in muon decay. We find that the correction can be described, within good 
approximation, by a hnear function in the electron energy. Explicit expressions for 
the form factors needed in an approach based on dispersion integrals are given. 
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At present, the best value of the muon hfe time is = (2.19703 ± 0.00004) fis [0. In 
fact, experimental data have reached such a precision that quantum corrections have to 
be taken into account when comparing this experimental result with theoretical predic- 
tions. To match this accuracy from the theory side, two-loop radiative corrections to the 
muon decay in the full electroweak Standard Model are needed. To perform the required 
calculations is a formidable, but not impossible, task. A step in this direction is the cal- 
culation of the purely electromagnetic corrections to the total decay rate at order in 
the Fermi theory, which has been performed by van Ritbergen and Stuart 0, ^ (see also 
0). Besides that, also the 0{Nfa'^) corrections in the Standard Model have been con- 
sidered in Refs. |[. The decay spectrum is known to order a since long; however, the 
corresponding corrections for polarized muon decay were calculated only recently in [0. 
The leading logarithmic QED corrections of order O (^a^ ln^(m^/me)j to the muon decay 
spectrum were considered in Ref. 

The spectrum calculation is different from the one for the life time since the Kinoshita- 
Lee-Nauenberg theorem |^ is not in effect. Consequently, powers of the large logarithm 
ln(m^/me) do not cancel in the calculation of the electromagnetic corrections to the 
spectrum. This becomes obvious when fitting the spectrum corrected to order 0{a) to 
the Michel spectrum: the resulting effective Michel parameter differs by about 6 % from 
its lowest-order value [10], a correction which is more than 10 times larger than the 
corresponding correction to the muon life time. At order 0{a^) the radiative corrections 
may be expected to be of the order of several per-mille and will become important for 
future high-precision experiments. 

Given this perspective we review the calculation of the hadronic contribution to the 
energy spectrum of the final-state electron in muon decay. This contribution is not ex- 
pected to be logarithmically enhanced, but nonetheless is required for an eventual com- 



plete second-order calculation. Further details of the calculation are given in Ref. ||TT 
We consider the decay of a muon in its rest system, 

fi~{p) -> e~{p') + u^{qi) + Ve{q2) , (1) 

and define momenta as shown in (|l]). It is convenient to introduce the dimensionless 
variable x = 2Ee/ m to denote the ratio of the energy of the decay electron with respect 
to the muon mass = m. Neglecting the electron mass, p'"^ = 0, the kinematically 
allowed range is < x < 1, and the momentum transferred from the charged particles to 
the neutrino pair, q = qi + q2 = p — p' is determined by = (1 — x)m? . 

The matrix element M. for ([l|) in the Fermi theory can be calculated most conveniently 
after a Fierz rearrangement factorizing the amplitude into a current which describes 
the /xe transition and a current for the v^Ve interaction. After squaring and summing 
(averaging) over spins, one can write as a product of two corresponding tensors. 

The one pertaining to the neutrino interaction can be integrated over the unobserved 
momenta of the neutrinos independently, leading to 

= q,,q,y - Qi^uQ^ ■ (2) 
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Figure 1: Feynman diagrams describing a self energy insertion in the photonic one-loop 
corrections to the fie vertex. 



This tensor will be contracted with C^i, = J*Ju, with = Ue(p')A^(g)M^(p), where 
A^(g) is the effective vertex of the four-fermion interaction. At the lowest order, A^(g) is 
identified with 



Aj = ^7.(l-75), 



(3) 



where Gp = (1.16637 ± 0.00001) x 10"^ GeV~^ is the Fermi coupling constant. 

Hadronic contributions to the radiative corrections to the current at order 0{a'^) are 
described by the Feynman diagrams shown in Fig. [^. The hadronic vacuum polarization 



P + iO 



(4) 



is inserted in a one-loop vertex correction. The vacuum polarization can be related to the 
measured cross section for e^e~ hadrons with the help of a dispersion relation 



nhad(^2^ = ^ r '^R^s)—-^ , R{S) 



0" s; e ' e 



hadrons) 



(5) 



where the integration starts at the two-pion threshold, Sthr = 4m^. Therefore, the calcu- 
lation corresponds to a one-loop vertex with a photon of mass y/s, i.e. using a propagator 



—1 



fc2 - s + iO 

The result can be written in the form 

a r°° ds 



kfxky 



A,(g) 



Svr Jsthi s 



(6) 
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where the vector function A„ can be decomposed into Lorentz-covariants as 



^ ^ujRg+{s\q ,m)^^ujR9-\s\q ,m) , (8) 
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with ujR^L = (1 ±75)/2. The calculation is straightforward and corresponds to that of the 
one-loop vertex correction in QED, with the difference that (i) the exchanged "photon" 
is massive with mass ^/s, (ii) the coupling is purely left-handed, and {Hi) the two fermion 
lines have different masses. In fact, except for small the electron mass can safely be 
neglected. After contraction with the neutrino tensor iV^,y, Eq. (Q), only the form factors 
/(s; g^, m^) and g^, m^) will remain in the final result. 



Explicit results for the relevant form factors / and are |TT 
f{s;q^,m^) = -^<! 2(m" - g" - s) 1 + 
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where Co = Co(m^, 0, g^; rn? , s, 0) is the three-point integral (cf. Fig. 0a), whereas Bq and 
A denote the tadpole and two-point integrals [|1^, respectively (see Eqs. (p^-(p^ 
below). The complete set of results, including that for g^, is presented in ||TT| . 

Self-energy diagrams (Fig. |I|b,c) contribute to the coefficient of 7^ only. The result is 



/sE(s;g^m^ 
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where the last term, |, comes from the self energy on the massless (electron) leg. Using 
recurrence relations [1^ (see also Appendix A of [0), the derivative in (|lT]) can be 
represented as 



dBo{p'^; 17?, s) 
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The required tadpole and two-point integrals are [0, ^ 
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where (3 = \Jl — Am? / and /xdr is the scale parameter of dimensional regularization. In 
Eqs. (|T2D-(|T4D, terms containing A = 1/e — Invr — 7e represent the ultraviolet singularities 
which cancel in the final results (PD-(p!0|). 

Finally, we need the three-point scalar function Cq for positive values of q^. For 
arbitrary values of and me, it can be presented in the following symmetric form: 
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where 



l+/3e 



Ami 



ml + m^2 — q^ — ^ X{ml, m^^, 



\ ml + m^2 - q2 jf. ^\{ml,m^^,q 



(16) 
(17) 



and \{x, y, z) ^ + + — 2xy — 2yz — 2zx is the standard notation for the Kallen 
function. This expression is regular in the limit m^ 0. A compact result for this case 
is given in Eq. (34) of |jTT| . 

Performing the dispersion integral one obtains the form factors 



a r°° ds 



a /■°° ds 



f{q\ m') = — —R{s)f{s- q\ m') , g+{q\ m') = — —R{s)~g+{s- q\ m') . 

(18) 

Since we are interested in the 0{a'^) correction, it is sufficient to keep only terms of first 
order in / and in the decay spectrum which then can be written in the form 



1 dF 



1 dF 



--^2x^[(3-2.)(l + 2/(.))+.,4^)]^-- 



(l + r(a;)). 



(19) 
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= 2x^(3 - 2x)ro , r(x) = 2/(x) + -^g^{x) , (20) 
„ 6 — 2x 

Born 



where /(x) = /((! — x)rn? , m^) and g+{x) = (7+ ((1 — x)m^, m^). 

For our purpose, the function R{s) describing the hadronic cross section of e+e^ anni- 
hilation can be modeled by a combination of experimental data and analytical results from 
perturbative QCD. Since we are going to calculate a small correction, it is not necessary 
to invoke the most sophisticated treatment as needed, for example, when calculating the 
hadronic contribution to the fine structure constant a{mz)- At low s < 2.5 GeV^ we use 
experimental data from ALEPH parametrized in |jl2| or provided directly by ALEPH [jl3| 



from a measurement of the isovector r spectral function. These data are complemented 
by the resonance contributions from the isospin-0 light mesons u and 0. Above s = 2.5 
GeV^ we use the QCD prediction for R{s) due to light quarks at order 0{as)- 

Since the data in the cc-channel published by various groups are in a large part of 
the energy range inconsistent, we apply in this case the QCD-based approach of analytic 
continuation by duality ^M. The data region can be chosen to extend only over the sub- 



threshold resonances, i.e. one can calculate the contribution coming from the cc-channel by 
a combination of data describing the J/\E'(1S') and J /'^{2S) resonances and the prediction 
of perturbative QCD. We checked that the results obtained this way are consistent with 
those of the standard approach using the new BES data |T5[] . 

The correction to the total decay rate, AF = dx (dF/dx), was calculated before in 
0. Our result, 

AFhad=^ -0.0421 (a/7r)'Fo, (21) 

agrees perfectly with the corresponding number —0.042 given in p[. The resulting correc- 
tions to the spectrum ([19|) are shown in Fig. |^. At small x, the corrections are positive, 
but the correction to the total decay width is dominated by the negative values at x^O.18. 
The dependence of the form factors on x is to a very good approximation linear: 



/(x) - (0.0071 - 0.0378x) (a/7r)^F 



; 



g+{x) - -0.0067 (a/7r)'Fo, (22) 



r(x) ~ (0.0148 - 0.0813x)(a/7r)^Fo. 

For 5f+, the coefficient of the term linear in x is very small and is therefore omitted. Note 
that this behaviour cannot be described by a simple redefinition of the Michel parameter. 
Since Gi? is a free parameter in the Fermi theory, the correction to the total decay width 
is not observable; it can be absorbed by a suitable redefinition of the Fermi constant. 
However, the modification of the spectrum is, in principle, observable. 

The correction to the total decay rate can be split up into the various contributions to 
the hadronic vacuum polarization, as shown in Table |1]. The form factor / of the 7^ term 
contributes —0.0387, whereas the correction due to (7+, —0.0033, is smaller by one order of 
magnitude. The total correction (the contributions due to / and g^) is saturated to 81 % 
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Figure 2: Results for the form factors f, and r defined in ^TB^), ([^ 



(80.8% and 89%, respectively) by the contributions from small s below 2.5 GeV^. Only 
5 % of AFhad is due to charmed states, and the bottom sector is completely negligible. 

The numerical results given here take into account the 0{as) QCD corrections in R{s). 
Using the leading-order expression for R{s), i.e. without the correction factor (l + as/vr) in 
the light-quark contribution at large s and the charm-quark contribution, the final result 
would be —0.04149, i.e. changed by 1.4%. We conclude that a more refined treatment 
which would include higher orders of perturbative QCD and mass-dependent corrections 
in the heavy-quark sector is not required for our purpose. 

The evaluation of the dispersion integrals has been performed using standard numerical 
integration routines up to a value Smax of several hundred GeV^. The contribution above 
this value was obtained with the help of Maple using the asymptotic expansion of the 
form factors for large s (see in the Appendix of Ref. |]TT|). For intermediate values of s, 
good consistency of both procedures has been verified. 

The same set of formulae can be used to calculate the contributions from a fi~^fi~ or 
a r+r" loop insertion. We find that the tau loop gives a very small contribution, about 
1.5 % of the one from the muon loop, in agreement with Ref. 0. Therefore we give only 
results for the muon loop, where one has to insert in Eq. ([T^) 




and Sthr 



4m^ 



(23) 
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Contributions 


to AFhad 




1 


< s < 0.2 GeV2 


-0.00129 


3.1 % 


2 


UJ 


-0.00223 


5.3 % 


3 




-0.00264 


6.3 % 


4 


0.2 < s < 2.5 GeV^ 


-0.02804 


66.6 % 


5 


s > 2.5 GeV^ 


-0.00564 


13.4 % 


6 


7/^(15) 


-0.00066 


1.6 % 


7 


J/^(25) 


-0.00017 


0.4 % 


8 


charm, s > 4m^ 


-0.00138 


3.3 % 


9 


bottom, s > Ami 


-0.00003 


0.1 % 




Sum 


-0.04207 


100 % 



Table 1: Contributions to the corrections of the total decay rate 



With this input we obtain 

Ar„~ -0.0364 (a/7r)' To, (24) 

which perfectly agrees with the exact result given in Ref . . The results of a linear fit of 
the form factors for the muon-loop insertion are: 

/muon(a:) ^ (0.0130 - 0.0414a;) (a/vrf To, 
^+,muon(a;) ^ (-0.0090 + 0.0005x)(a/7r)^ To, (25) 
r^y^onix) ^ (0.0267 - 0.0898x)(a/7r)^ To. 

When analysing the similar contribution with an e~^e~ loop insertion (which represents 
one of the pure QED contributions), one should carefully treat logarithmic singularities 
occurring in the limit rrie — > 0. For rUe ^ 0, the general result (0) can be used. 

One can see that the energy spectrum in the decay of an unpolarized muon is corrected 
by a smooth function in x due to hadronic contributions at order 0{a'^). At both ends of 
the spectrum no particularly outstanding enhancement or suppression is observed. 

The calculations described in this paper constitute only the most straightforward part 
of a full calculation which would be necessary before the expected future high-precision 
data can be confronted with theoretical predictions. This will not only be necessary for a 
meaningful test of the electroweak Standard Model, but also when searching for physics 



beyond the Standard Model [16 
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